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ABSTRACT 

Super star clusters — extremely massive clusters found predominately in starburst environments — 
are essential building blocks in the formation of galaxies and thought to dominate star formation in 
the high-redshift universe. However, the transformation from molecular gas into these ultra-compact 
star clusters is not well understood. To study this process, we used the Submillimeter Array and 
the Plateau de Bure Interferometer to obtain high angular resolution (~ l'/5 or 160 pc) images of the 
Antennae overlap region in CO(2-l) to search for the molecular progenitors of the super star clusters. 
We resolve the molecular gas distribution into a large number of clouds, extending the differential 
cloud mass function down to a 5a completeness limit of 3.8 x 10 5 M Q . We identify a distinct break in 
the mass function around log M mo i/M Q ss 6.5, which separates the molecular clouds into two distinct 
populations. The smaller, less massive clouds reside in more quiescent areas in the region, while the 
larger, more massive clouds cluster around regions of intense star formation. A broken power-law fit to 
the mass function yields slopes of a = — 1.39± 0.10 and a = — 1.44 ± 0.14 for the low- and high-mass 
cloud population, well-matched to the mass function found for super star clusters in the Antennae 
galaxies. We find large velocity gradients and velocity dispersions at the locations of intense star 
formation, suggestive of compressive shocks. It is likely that these environmental factors contribute 
to the formation of the observed massive molecular clouds and super star clusters in the Antennae 
galaxies. 



Subject headings: galaxies: individual (NGC 4038/4039) 
galaxies: star clusters 



galaxies: ISM — galaxies: starburst 



1. INTRODUCTION 

In areas of lower extinction between the dust lanes of 
starbursts, observations by the Hubble Space Telescope 
(HST) have revealed an important connection between 
the starburst phenomenon and a class of objects known 
as "super star clusters" (SSCs; Larsen 2010; Portegies 
Zwart et al. 2010). SSCs represent the most extreme 
mode of star formation known. Aside from their extraor- 
dinary luminosities, which range from 1-100 times that 
of the R136 cluster of 30 Doradus in the Large Mag- 
ellanic Cloud (Holtzman et al. 1992; Whitmore 2003; 
Larsen & Richtler 2004), these clusters are exceptionally 
compact (radii < 1 — 5 pc; Larsen 2010; Portegies Zwart 
et al. 2010) and young (< few - 50 Myr; Portegies Zwart 
et al. 2010). Estimates based on high- resolution spec- 
troscopy and population synthesis models suggest that 
cluster masses range from 10 4 to 10 6 M (e.g., Ho & Fil- 
ippenko 1996a,b; Zhang & Fall 1999; Mengel et al. 2002; 
Larsen & Richtler 2004; Melo et al. 2005; McCrady & 
Graham 2007). SSCs comprise a significant fraction of 
the UV, Ha, and B-band light emitted by young stars in 
starburst regions (Barth et al. 1995; Meurer et al. 1995; 
Maoz et al. 1996; Zepf et al. 1999; Fall et al. 2005). This 
suggests that in these regions SSCs are the basic building 
blocks of star formation (Ho 1997). 

The similarity in size and mass of the SSCs and glob- 
ular clusters suggests that the SSCs might be present- 
day analogs of young globular clusters (Holtzman et al. 
1992; McCrady & Graham 2007). Previously thought 
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to form strictly in the primordial epoch of galaxy evo- 
lution, globular clusters may in fact be forming in the 
current epoch as SSCs in starburst environments (Ash- 
man & Zepf 2001). 

The formation of SSCs, however, is not well under- 
stood. Some key questions include: How is the mass 
of the precursor molecular cloud, equivalent to a giant 
molecular cloud (GMC), converted to stars quickly and 
efficiently enough for the star cluster to remain bound? 
How are the conditions in starbursts different from "or- 
dinary" star formation? Are there other distinguishing 
factors that differentiate progenitor clouds of SSCs from 
normal GMCs? 

Recent studies suggest that it is high environmental 
pressure that increases the gas density and ultimately 
leads to the formation of massive star clusters (e.g., Ash- 
man & Zepf 2001; Herrera et al. 2011). One hypothe- 
sis suggests that a sudden galactic-scale increase in ex- 
ternal pressure, possibly from shock compression associ- 
ated with mergers, crushes entire GMCs rapidly and ef- 
ficiently, forming individual SSCs (Jog & Solomon 1992; 
Ashman & Zepf 2001). CO (2-1) line observations of M82 
with the Owens Valley Millimeter Array (OVRO) with 
17 pc spatial resolution by Keto et al. (2005) support this, 
showing that the mass spectrum of GMCs in the star- 
burst region is indistinguishable from the upper range of 
the mass spectrum of SSCs in M82 (Melo et al. 2005). 
These observations further show that the clouds with ac- 
tive star formation are in supersonic compression, with 
both lateral and line-of-sight velocity gradients within 
the individual clouds. 

In this hypothesis, the compression itself is created by 
the rapid ionization of the interstellar medium (ISM) 



surrounding the progenitor GMCs. The initial ioniza- 
tion may be caused by the dissipation of kinetic energy 
in the diffuse ISM of colliding galaxies. The burst of 
massive star formation that follows then releases enough 
energy to further propagate the ionization and compress 
nearby molecular clouds. This hypothesis for ionization- 
driven implosion draws a fascinating connection between 
the present-day SSCs and the ancient globular clusters in 
that the globular clusters may also have formed in clouds 
compressed by rapid ionization in the epoch of reioniza- 
tion in the primordial universe at z k, 6 (Cen 2001; van 
den Bergh 2001). 

Alternatively, SSCs may form in the high-pressure en- 
vironment of the centers of super-giant molecular com- 
plexes (M > 10 7 M Q ; Harris & Pudritz 1994; Wilson 
et al. 2003). In this hypothesis, cores that are the size 
of GMCs, but are within larger molecular complexes, are 
slowly compressed (on a crossing time) by the pressure of 
the overlying molecular gas. Cluster formation proceeds 
as a scaled-up version of normal star formation in our 
own Galaxy (Elmegreen & Efremov 1997). In support 
of this alternative, CO(1-0) line observations of the An- 
tennae galaxies show huge concentrations of molecular 
gas that could be interpreted as super-giant molecular 
complexes (Wilson et al. 2003). Yet because the angu- 
lar resolution of these observations (3'.'2x4'.'9; ^420pc 
at the adopted distance of the 22 Mpc) is insufficient to 
resolve individual GMCs, the question remains whether 
the appearance of the extraordinary massive clouds is 
simply the result of limited spatial resolution. 

In this paper we revisit the Antennae galaxies, 
NGC 4038 and NGC 4039, with higher resolution imag- 
ing of the molecular gas to study the formation of SSCs. 
As one of the most spectacular starbursts in the local uni- 
verse, the Antennae has been studied in great detail over 
the last several decades in various wavelengths, ranging 
from the X-rays to the mm/cm regime, and everything in 
between (e.g., Mirabel et al. 1998; Neff & Ulvestad 2000; 
Wilson et al. 2000, 2003; Fabbiano et al. 2001, 2003; Hi- 
bbard et al. 2001, 2005; Zhang et al. 2001; Kilgard et al. 
2002; Mengel et al. 2002; Zezas et al. 2002; Zezas & Fab- 
biano 2002; Fall et al. 2005). It is also one of the best- 
studied SSC populations (Whitmore & Schwcizer 1995; 
Whitmore et al. 1999, 2010). 

The overlap region of the Antennae, where NGC 4038 
and NGC 4039 collide, is one of the most intensely 
and violently star-forming regions in the nearby universe 
(Evans et al. 1997; Mirabel et al. 1998; Gao et al. 2001; 
Gao 2008; Brandl et al. 2009; Klaas et al. 2010; Zhang 
et al. 2010). The star formation rate (SFR) for the entire 
NGC 4038/39 system is high, estimated to be between 
6.6-20 Moyr^ 1 (Zhang et al. 2001; Brandl et al. 2009). 
However, most of the star formation occurs in the overlap 
region, with SFRs of 0.7-2.0 il^yr" 1 for individual star- 
forming regions on size scales of 500 x 1200 pc 2 (Brandl 
et al. 2009). Recent studies also find intriguing evidence 
for shock-driven star formation (Gao et al. 2001; Zhang 
et al. 2010), in agreement with the Jog & Solomon model. 
Detailed population synthesis analysis of the clusters in 
these regions finds the median age of the clusters to be 
some of the youngest in the system, <10Myr in many 
cases (e.g., Snijders et al. 2007; Zhang et al. 2010). 

In this paper, Section 1 comprises this introduction. In 




Fig. 1. — : Three-color image of the Antennae from the HST 
Advanced Camera of Surveys (ACS), using data taken by 
Whitmore et al. (2010). Filters used to create the image are 
F435W (blue), F550M (green), and F814W (red). The image 
has dimensions of 160" x 195", and it is displayed using an 
arcsinh stretch (Lupton et al. 2004). The white circle shows 
the primary beam of our SMA observations, and red circles 
show our PdBI pointings. See Table 1 for more information. 

Section 2 we present new, high-resolution CO (2-1) data 
from the Submillimeter Array (SMA; Ho et al. 2004) and 
the IRAM Plateau de Bure Interferometer 3 (PdBI). Sec- 
tion 3 presents some basic properties of the molecular gas 
in the observed region and discusses the implications. In 
Section 4 we compare the distribution of molecular gas 
mass, velocity, and velocity dispersion with star forma- 
tion tracers. In Section 5, we summarize our results. 
Throughout this paper we adopt a distance of 22 Mpc 
for the Antennae (Schweizer et al. 2008). 

2. OBSERVATIONS 

2.1. The Submillimeter Array 

The SMA observations were taken in May 2008 and 
February 2011 in the compact and extended array con- 
figurations, respectively, with a single pointing with a 
phase center (a,5) J2 noo = (12 h 01 m 54 s .70, -18°53'05"0). 
Figure 1 shows our spatial coverage of the overlap re- 
gion 4 . At 230 GHz, the primary beam of the SMA covers 
~54". Baselines ranged from 16 to 139 meters (compact 
array) and 44 and 226 meters (extended array). Typi- 
cal system temperatures ranged between 100 and 300 K 
for both observations. The correlator was tuned to the 
frequency of CO(2-l) at 230.53799 GHz, with Doppler 

3 IRAM is supported by INSU/CNRS (France), MPG (Ger- 
many) and IGN (Spain). 

4 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space Tele- 
scope Science Institute. STScI is operated by the Association of 
Universities for Research in Astronomy, Inc. under NASA contract 
NAS 5-26555. 



tracking accounting for the recessional velocity of the 
Antennae galaxies. For the May 2008 track, the correla- 
tor was configured to cover a 2 GHz bandwidth with 24 
overlapping windows each 104 MHz in width, with a spec- 
tral resolution of 0.4 MHz (~0.5kms _1 channels). The 
February 2011 track was configured slightly differently 
due to the bandwidth doubling upgrade by the SMA in 
2009, covering 4 GHz with 48 overlapping windows each 
104 MHz in width, with a spectral resolution of 0.8 MHz 
(^lkms -1 channels). 

The SMA data were reduced using the MIR data re- 
duction package developed at OVRO and the SMA. Data 
were flagged for bad channels, antennas, weather, and 
pointing. Bandpass calibration was done on a bright 
quasar, 3C 454.3 for the May 2008 observation and 3C 84 
for the February 2011 observation. We performed atmo- 
spheric phase calibrations with 1127—189 and 1256—057 
every 10 minutes. Absolute flux calibration was done 
using Callisto and Titan. The channels were Hanning 
smoothed and averaged to achieve a final resolution of 
4.9 km 



,-i 



2.2. Plateau de Bure Interferometer 

The PdBI observations were taken between February 
and April of 2006 in the B and C arrays, respectively. 
Because of the smaller (~21") primary beam of the 
PdBI 15 meter antennas, three separate pointings were 
obtained to cover approximately the same area of the 
overlap region as the SMA primary beam at 230 GHz 
(Figure 1 & Table 1). Each pointing was observed with 
a single ^7 hour track at the PdBI. Baselines ranged 
from 24 to 82 meters in the C array and 40 to 280 meters 
in the B array. Typical system temperatures ranged from 
300 to 1000 K due to the low elevation of the Antennae 
for the PdBI. The correlator was configured to cover the 
CO(2-l) line with four windows each 160 MHz in width, 
with 160 channels in each window with spectral resolu- 
tion of 1.25 MHz (~1.6kms _1 channels). 

The PdBI data were reduced using the GILDAS CLIC 
package. Bandpass and absolute flux calibration were 
done using 3C 273, and the atmospheric phase was cali- 
brated using 1124—018 every 20 minutes. The channels 
were Hanning smoothed and averaged to achieve a final 
resolution of 4.9kms~ 1 . 

2.3. Imaging 

To recover the more extended flux from the shorter 
spacings and increase the signal-to-noise ratio, we com- 
bine the high-resolution SMA and PdBI data together 
with the SMA compact array observations in the (u,v) 
plane. After reduction, the PdBI (u,v) data were saved 
as FITS files, then converted into the MIRIAD format 
using the task fits. Header information including tele- 
scope name, longitude, latitude, and channel resolution 
were manually put into the PdBI MIRIAD datasets. We 
combined the data with the MIRIAD software package 
(Sault et al. 1995), using natural weighting and an ad- 
ditional weighting in inverse proportion to the noise as 
estimated by the system temperature. The data cube 
was cleaned to a cutoff of 2 er c han in the residual image. 
The synthesized beam of the combined SMA+PdBI im- 
age is 3'.'3xl'.'5, with a position angle of 41°. The rms in 
a single 4.9kms _1 channel (er c han) is 30mJybeam _1 . 





TABLE 1 


Telesc 


ope Pointing Centers 


Position 


Phase Center 




("2000 , <52000 ) 


SMA Comp 


12 h 01 m 54 s .70,-18°53'05'.'0 


SMA Ext 


12 h 01 m 54 s .70,-18°53'05'.'0 


PdBI A 


12 h 01 m 55 s .00,-18°52 / 50 / /0 


PdBI B 


12 h 01 m 54 s .70,-18°53'05 / . / 


PdBI C 


12 h 01 m 54 s .00,-18°53'15y0 



Velocity- integrated emission (zeroth moment), 
intensity- weighted velocity (first moment), and intensity- 
weighted dispersion (second moment) maps were made 
from a masked data cube. The masks were made by 
smoothing the original cube to a resolution of 6" and 
masking emission that was less than la or did not 
cover two contiguous channels. Figure 2 shows the 
velocity-integrated CO(2-l) emission of our combined 
PdBI+SMA observations of the overlap region. 

3. PROPERTIES OF THE MOLECULAR GAS 
3.1. Finding Clouds 

To extract individual clouds from our data cubes, we 
utilize the algorithm clumpf ind by Williams et al. (1994) 
on our data cube where the PdBI and SMA primary 
beams overlap, clumpf ind contours the input data cube 
in steps that are integer multiples of the input channel 
rms and identifies clouds as emission peaks separated by 
more than one contour. The extent of each cloud is de- 
termined by following each peak down to the first (low- 
est) contour. Multi-peaked emission is separated by a 
"friends-of-friends" algorithm where pixels adjacent to 
the peak and subsequent adjacent pixels are identified as 
related. 

The default contour and detection thresholds recom- 
mended for clumpfind by Williams et al. (1994) are 
2cr c han- We tested three different detection thresholds 
(2, 2.5, and 3er c han while holding the contour threshold 
constant at 2<7 c h a n)- However, we find that most of the 
results do not depend strongly on the choice of detection 
threshold. We discuss this in more detail with the indi- 
vidual results. We set the clumpfind search parameters 
to contour at twice the channel rms (2 <r c han) with the de- 
tection threshold at 3 CT c han to detect the faintest clouds 
with a higher confidence level and better avoid spurious 
detections of low-mass clouds. 

3.2. Fluxes and Mass Estimates 

We estimate the mass of each cloud from its CO lumi- 
nosity following Solomon et al. (1992): 
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where F. qq is the velocity-integrated flux, ^ rost is the rest 
frequency of the line in GHz [230.538GHz for CO(2- 
1)], z is the redshift, and Dl is the luminosity distance, 
which we fix to 22Mpc (Schweizer et al. 2008). The CO 
luminosity is related to the molecular gas mass by 
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Fig. 2. — : Velocity-integrated CO(2-l) emission from the combined SMA+PdBI data of the overlap region of the Antennae. 
Contours start at 3 times the rms of the image (a — 3.71 Jybeam -1 kms -1 ) and increase by 3a. Intensities range from 0-100 
Jybeam -1 kms -1 . The synthesized beam of this image is 3'.'3xl'.'5, corresponding to a physical scale of 374x 160 pc. We denote 
the locations of the clouds found by clumpf ind with circles. Symbol size, thickness, and color reflect the molecular gas mass of 
the cloud estimated from the CO brightness. We find multiple overlapping clouds, with the most massive clouds concentrated 
at the peaks of integrated CO emission and less massive clouds in the surrounding regions. Numbers indicate the clouds whose 
intensity-weighted velocity is shown in Figure 5. We discuss this figure in further detail in § 3. 



where aco = 2toh^co 

There has been some debate over the value of Xqo 
in the Antennae. Zhu et al. (2003) find Xqo to be 5- 
13 times smaller than the Galactic value (ATco,gai = 
2.8xl0 20 (Kkms -1 ) -1 ). Analysis of resolved CO (1-0) 
clouds by Wilson et al. (2003), however, suggest that 
Xqo i n the Antennae is much larger, with an average 
value of 1.3xl0 20 (Kkms -1 ) -1 and within the errors of 
the Galactic value. More recently, models of CO and 
[C II] data by Schulz et al. (2007) suggest that X co 
does not deviate much from the Galactic value. Based 
on these results and for consistency with Wilson et al. 
(2003), we adopt a CO =4.8M Q (Kkms -1 pc 2 ) -1 , which 
corresponds to a CO-to-Hi2 conversion factor Xqo = 
3xl0 20 (Kkms -1 ) -1 . 

We also assume a line intensity ratio of 
JcO(2-i)Aco(i-0) = 1-0. Zhu et al. (2003) find 
the /co(2-i)/^co(i-o) ime intensity ratio to range 
from 1.0 to 1.2 in the overlap region of the Antennae 
with their single-dish observations. This is also the 
typical value found for nearby disk galaxies, as well 
as the central region of the Milky Way (Braine et al. 



1993; Sawada et al. 2001). We multiply M H2 by 1.36 
to account for heavier elements to obtain the total 
molecular gas mass, M mo \. 

One concern, especially with the higher resolution 
imaging, is the problem of missing flux from the lack 
of short spacings in interfcrometric observations. In or- 
der to estimate the percentage of the total flux captured 
by our interferometric observations, we first estimate the 
total flux in the region covered by other observations. 
CO(1-0) single-dish mapping of the Antennae with the 
NRAO 12 m by Gao et al. (2001) gives a total flux of 
3172Jykms -1 . Note that the region mapped by Gao 
et al. (2001) is much larger than the region covered by our 
observations or that of Wilson et al. (2003). They used 
the single-dish spectra of Gao et al. (2001) to estimate 
a total CO(1-0) flux of 1654 Jy kms -1 in their smaller 
area mapped by the OVRO interferometers. We can use 
the estimate of Wilson et al. (2003) to determine the 
total CO (1-0) flux in the area mapped by our observa- 
tions. Since our observations do not include NGC 4038, 
we subtract 256 Jy kms -1 from Wilson's estimate to ob- 
tain a total single-dish CO(1-0) flux for our region of 



lSgSJykms- 1 . 

We can compare the predicted total CO flux with the 
flux measured in our interferometric observations. We 
determine the flux by summing up the fluxes of the clouds 
found by clumpf ind in our combined SMA+PdBI data. 
We obtain a total flux of 1358Jykms _1 , which is very 
close to the single-dish flux assuming a line intensity ratio 
of Jco(2-i)/^co(i-o) = 1-0. If we assume a line intensity 
ratio of 7co(2-i)/^co(i-o) = 1-1 (Zhu et al. 2003), then 
we are missing only ~12% of the total flux. This indi- 
cates that most of the CO mass may be in discrete clouds 
rather than distributed in the diffuse ISM in the overlap 
region. Note that this may not be the case for regular 
star-forming disks, where the gas tends to be more dif- 
fuse. Lower cutoff thresholds in clumpf ind yield total 
fluxes of 1472 and 1614Jykms -1 for the 2.5 and 2cr c h a n 
thresholds, respectively, which imply that noise clumps 
may be picked up as real clouds at these threshold levels. 

3.3. Basic Cloud Properties: Two Populations 

With the clumpf ind parameters adopted in § 3.1, we 
detect a total of 132 clouds. We list all the clouds and 
their basic properties in Table 2. As expected, if we 
use the lower detection thresholds, clumpf ind extracts 
a larger number of clouds (total of 332 and 224 clouds 
for the 2 and 2.5<7 c h a n thresholds, respectively). Most of 
the new clouds detected with the lower thresholds are 
low mass and unresolved — the number of clouds with 
log M mo i/M© < 6.5 doubles and triples with the 2.5 and 
2c c han thresholds and therefore are not significant in our 
analysis. 

We find a large range of values for the basic prop- 
erties of the 132 clouds (Figure 3 and Table 3). The 
molecular masses (discussed in § 3.2) span over 3 or- 
ders of magnitude in dynamic range, from M mo \ w 10 
to 2.4 x 1O 8 M , with a mean of 2.0 x 1O 7 M , a me- 
dian of 3.0 x 1O 6 M , and standard deviation (calcu- 
lated in log space) of 1.0 dex. One of the most inter- 
esting characteristics of the mass distribution is that it 
appears distinctly bimodal, with a maximum around our 
sensitivity limit of M mo \ « 10 5 M and clear secondary 
maximum near M mo \ « 1O 7 M . The exact location of 
the minimum is somewhat arbitrary, within the range 
of log M mo i/M Q w 6.3 — 6.7. From considerations of 
the mass function (§ 3.7), we pick a fiducial value of 
logM mo i/M = 6.5 to separate the two populations. 
With this choice, there are 68 clouds in the low-mass 
group and 64 clouds in the high-mass group, whose re- 
spective median masses are M mo \ — 2.3 x 10 5 M Q and 
M mo i = 1.4xlO 7 M . 

Values of radius and velocity dispersion estimated by 
clumpf ind are corrected for spatial and spectral resolu- 
tion following equation A7 from Williams ct al. (1994). 
We find that the radii of the clouds span a very large 
range, from r — 35 to 317 pc. Of the 132 clouds found, 
64 are unresolved spatially. We assign an upper limit 
of 80 pc to the radius of these clouds, half of our mini- 
mum resolution diameter of 160 pc. Using the Kaplan- 
Meier product-limit estimator (Feigelson & Nelson 1985), 
which properly accounts for the upper limits, we find 
a mean radius of 119.2±7.3pc and a median value of 
77.0 pc for the entire sample. The distribution of sizes, 
like that of mass, appears bimodal. In fact, given the 



large fraction of upper limits at the small-r end of the 
distribution, the true separation between the two peaks 
must be even greater than it appears. Splitting the 
sample into the low-mass and high-mass groups defined 
above, the respective mean and median radii become r = 
(57.0±7.1, 48.0) pc and (180.6±8.8, 170.0) pc. In light 
of the large number of upper limits, the statistics for the 
low-mass subsample should be treated with caution. 

The clouds on the upper end of the size distribution are 
larger than any known GMCs in the Milky Way (Solomon 
et al. 1987; Heyer et al. 2009) or in nearby disk galaxies 
(Bolatto ct al. 2008). While blending may be an issue 
(but see arguments below as to why we do not think this 
is a major problem), at face value the largest clouds rep- 
resent a distinct class. Wilson et al. (2000) described 
these very large clouds as super-giant molecular com- 
plexes. Our spatial resolution of 80 pc is high enough to 
resolve the most massive clouds into smaller clouds, but 
nonetheless they appear to be single entities. There is 
very good correspondence between cloud size and cloud 
mass; the most massive clouds are also the largest in 
physical extent. The corollary implication, therefore, is 
that the most massive clouds in our sample are intrin- 
sically as massive as they appear. Their exceptionally 
high masses are not a spurious artifact of insufficient res- 
olution. This is one of the central results of this study. 

Of the 132 clouds found, 86 are resolved in veloc- 
ity, with velocity dispersions (o~ v ) ranging from 0.8- 
25.6kms~ 1 . For the unresolved clouds, we assign an 
upper limit of 4.9 km s -1 , the channel size of our data 
cube. There is an overlap of 45 clouds that have limits 
in both radius and velocity dispersion. Note that a few 
clouds have radius and/or o~ v smaller than the upper limit 
because they were bright enough to be resolved. Unlike 
the sizes and masses, the distribution of velocity disper- 
sions is not as cleanly separated between the high-mass 
and low-mass groups, although lower mass clouds tend to 
have lower values of o~ v . We calculate a mean and median 
velocity dispersion of 7.8±0.6kms _1 and 6.2 kms" 1 , re- 
spectively, for the sample as a whole. Considered sep- 
arately, the low-mass and high-mass groups have mean 
and median dispersions of a v = (2.6±0.3, 2.0) kms^ 1 and 
(13.3±0.6, 12.0) kms -1 , respectively. Again, the statis- 
tics of the low-mass group may be uncertain because of 
the large fraction (68%) of upper limits. 

3.4. Molecular Gas and Cloud Kinematics 

We show the intensity-weighted velocity and velocity 
dispersion map of our combined SMA+PdBI CO(2-l) 
data in Figure 4. Consistent with previous interfero- 
metric studies of molecular gas in the Antennae (Gao 
ct al. 2001; Wilson et al. 2003; Ueda et al. 2012), we find 
that the kinematics of the molecular gas in the over- 
lap region are extremely complicated. Several of the 
CO(2-l) intensity peaks, especially NGC 4039 in the 
lower right corner, show fairly obvious velocity gradients. 
Intensity-weighted velocity dispersions are large — often 
>50kms~ 1 — most likely due to multiple clouds along 
the line of sight, which we will discuss later. This does 
not contradict our conclusion that we can resolve individ- 
ual clouds because clumpf ind differentiates by velocity 
as well as spatially. 

The velocity map gives the average velocity at any 
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Fig. 3. — : Histogram of radius, velocity dispersion, and molecular gas mass for the clouds found by clumpf ind. Clouds that 
were unresolved spatially or in velocity are assigned upper limits of 80 pc and 4.9 km s -1 , respectively, which we mark in the 
first two panels with dashed vertical lines. The clouds are broken up into two groups, above and below log M mo i/M0 = 6.5 
(shaded and open histogram, respectively), which corresponds to a distinct break in the differential mass function discussed in 
5 3.7. 



TABLE 3 
Cloud Statistics 



r ff v log M mol 

Mean Error Median Mean Error Median Mean StDev Median 

(pc) (bus- 1 ) (log Me) 



All Clouds 

log M mol /M Q > 6.5 

log M mol /M Q < 6.5 


119.2 
180.6 

57.0 


7.3 

8.8 
7.1 


77.0 
170.0 

48.0 


7.8 

13.3 

2.6 


0.6 

0.6 
0.3 


6.2 

12.0 
2.0 


7.3 
7.6 
5.7 


1.0 

0.5 
0.4 


6.5 

7.2 
5.4 


-1.44±0.14 
-1.39±0.10 


4.01±1.01 
3.02±0.59 



Note. — Basic statistical properties (radius, velocity dispersion, and molecular gas mass) for the clouds extracted by clumpf ind. 
We break up the cloud statistics into two groups, log M mo i/M0 > 6.5 and log M mo i/MQ < 6.5, to show that there are significant 
differences between the general properties for these two populations. For r and a v , we list the mean, the error on the mean, and the 
median as calculated using the Kaplan-Meier product-limit estimator to account for the upper limits. The last two columns show the 
results from the ordinary least squares fits to the differential mass function for the two populations (Figure 9) in log-log space, where 
log d7V/dM mol =/3 + alog M mol . 



given position, weighted by the CO intensity, so veloc- 
ity information could be blurred out if there are multiple 
clouds along the line of sight. To preserve some velocity 
information, we made individual first moment maps for 
the 50 brightest clouds extracted by clumpf ind, using 
the masks provided by the algorithm. Figure 5 shows 
the first moment velocity maps and spectra for the 10 
most massive clouds found by clumpf ind in our com- 
bined SMA+PdBI data, with their spatial locations la- 
beled in Figure 2. We find that about half of the clouds 
show a velocity gradient indicating shear or rotation. 
Cloud 10 appears to have multiple velocity components. 
However, because the velocity peaks are separated by less 
than 2<7 c han, clumpf ind assigns the emission to a single 
cloud. 

3.5. Evidence for Gravitationally Bound Clouds 

We estimate the virial mass for each cloud from the 
resolution-corrected radius and velocity dispersion. The 
virial mass is given by 



M v 



5r< 
a g G 



(3) 



where a g is a geometric factor that depends on the den- 
sity profile. For a self-gravitating cloud in virial equilib- 
rium, we expect the density profile to scale as r~ 2 (Bo- 



denheimer & Sweigart 1968; Shu 1977; Keto & Caselli 
2010) and a g — 5/3. We plot the relationship between 
molecular gas mass and virial mass in Figure 6. The 
detections alone give (M v i r /M mo \) = 1.2, but this ratio 
drops to (Mvir/Mmoi) = 1.03 if we properly take the up- 
per limits into account. 

As mentioned in § 3.2, there is some debate over 
whether Xqo is much smaller in the Antennae (Zhu et al. 
2003; Wilson et al. 2003; Schulz et al. 2007). We show in 
the upper left corner of Figure 6 the distance each point 
would move to the left if Xqo was 2 and 5 times smaller, 
and to the right if Xqo was 2 times larger than our as- 
sumed value of Xco = 3xl0 20 (Kkms -1 ) -1 . Consider- 
ing the many uncertainties inherent in the determination 
of either M v ; r or M mo i, the data are entirely consistent 
with M v ; r /M mo i « 1. In other words, the clouds are 
gravitationally bound structures. This is a key point. It 
implies that the extreme properties of the largest clouds 
— their large sizes, internal velocities, and masses — are 
unlikely to be artifacts of spatial or velocity blending. If 
they were, then it would be a coincidence that the radii 
and velocity dispersions of the clouds would conspire to 
have just the right combination to produce M v ; r w M mo \. 

If we run clumpf ind with lower thresholds, we still find 
a one-to-one correspondence between molecular gas and 
virial mass, but some of the low-mass clouds fall quite a 
bit below this relation. These clouds typically have a v < 
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Fig. 4. — : Intensity-weighted velocity (left) and velocity dispersion (right) maps from our combined SMA+PdBI CO(2-I) 
data, masked at the 3 a level of CO intensity. White contours demarcate velocities and velocity dispersions in steps of 20 and 
15kms _1 , respectively. 



2 kms -1 and upper limits on radius. There are only one 
or two clouds found using the 3cr c han threshold with the 
same behavior, which reinforces our earlier conclusion 
that some of the additional low-mass clouds extracted 
using the 2 and 2.5er c han detection thresholds may not 
be significant above the noise. 

3.6. Cloud Scaling Relations 

Larson (1981) described the relationship between the 
size, linewidth (velocity dispersion), mass, and density of 
molecular clouds. Here we examine whether our clouds 
fall on these scaling relations and how they compare to 
typical GMCs from non-starburst galaxies. 

We plot the relationship between radius (r) and veloc- 
ity dispersion (a v ) for our clouds in Figure 7. Solomon 
et al. (1987) found a relationship of the form a v ex r ' 5 
for Milky Way GMCs. This correlation has been shown 
to hold even down to very small scales within individ- 
ual GMCs (Heyer et al. 2009). We find that most of 
our clouds fall on the r-a v relation, extending it to 
larger sizes and higher velocity dispersions. However, the 
largest clouds in our study appear to have systematically 
higher velocity dispersions than predicted for Milky Way 
GMCs. The clouds in the Antennae with r > 100 pc 
have velocity dispersions higher by a factor of ~ 2 com- 
pared to Galactic GMCs. Clouds found by Ueda et al. 
(2012) in their CO(3-2) observations of the Antennae, 
which were also obtained with the SMA and at a resolu- 
tion comparable to ours, also lie above the Galactic r— o~ v 
relation. In fact, the CO(3-2) dispersions are even larger 
than those measured in CO (2-1). For the 35 clouds in 
common between our sample and that of Ueda et al., 
(a v (3-2)/a v (2-l)) « 2. Matching our CO(2-l) clouds 
with the CO (1-0) clouds from Wilson et al. (2003) is 
difficult given the difference in resolution (3'.'2x4'.'9 com- 
pared to 3'.'3xl'.'5), but we note that most of the CO(l- 
0) clouds have velocity dispersions <10kms _1 , signifi- 
cantly smaller than the velocity dispersion observed in 
the CO(3-2) and CO(2-l) clouds. 

The large offset in the Ueda et al. velocity dispersions 



are especially curious. The typical size-linewidth relation 
in Galactic (Larson 1981; Solomon et al. 1987) and extra- 
galactic (Bolatto et al. 2008) clouds predicts that clouds 
found by higher density tracers such as CO (3-2) should 
have smaller velocity dispersions if they trace smaller size 
scales, precisely the opposite of what Ueda et al. (2012) 
observe. 

There are several possible explanations for the offset in 
velocity dispersion for the largest clouds. Assuming that 
the Milky Way GMCs represented by the solid line in Fig- 



M v 



ex ra,. 



( '- 5 ) found by Solomon 
that all the Milky 



ure 7 are in virial equilibrium (M mo ) 
the size-linewidth relation (a v ex r 1 
et al. (1987) implies that M oc r 2 
Way clouds have similar gas surface densities. For the 
GMCs in the Solomon et al. sample, this corresponds to 
a gas surface density E = 170 M Q pc~ 2 . If the Antennae 
clouds are in virial equilibrium, as they appear to be (Fig- 
ure 6), then one possibility is that they have higher sur- 
face densities than Milky Way GMCs. Gao et al. (2001) 
find that the overlap region has exceptionally high HCN 
(a high gas density tracer, > 10 4 cm -3 ) emission, equiv- 
alent to half of the total HCN emission in the system. 
Furthermore, the overlap region is extremely bright in 
the far- infrared (FIR; Klaas et al. 2010), which is known 
to be well-correlated with HCN (Gao & Solomon 2004). 
This correlation has been shown to hold for a wide range 
of size scales, from dense Galactic cores to extreme, high- 
redshift starbursts (Gao & Solomon 2004; Wu et al. 2005; 
Gao et al. 2007). 

We can roughly estimate the average H2 number den- 
sity in each cloud from its CO brightness, assuming the 
mass is distributed uniformly within a sphere. Larson 
(1981) found that the average number density ((n(H2))) 
varies inversely as the cloud diameter L [(n(H2))(cm 3 ) oc 
L(pc) -11 ], implying that the column density is essen- 
tially independent of size. Figure 8 shows that our clouds 
follow this relationship. At a given size, however, our 
clouds appear to be offset toward higher number densi- 
ties (by ~ 0.3 dex) than the molecular clouds from Lar- 
son (1981), as well as GMCs from Bolatto et al. (2008). 
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Fig. 5. — : Intensity-weighted velocity maps (middle panels) of the 10 most massive clouds extracted by clumpf ind (labeled 
in Figure 2), showing velocity gradients indicative of shear or rotation in many of the clouds. Clouds are numbered according 
to decreasing CO intensity. The angular size of each panel is the same (8.3" x 8.3"), but note that the color-scale is different 
for each panel, chosen to best show the velocity gradient. The white contours demarcate velocities in steps of lOkms -1 for all 
panels. The CO spectrum of each cloud is plotted in the top/bottom panels, with channel widths of 4.9 km s -1 . 



This is consistent with the argument for higher surface 
densities, and could also be the explanation for the even 
larger offset of the Ueda et al. (2012) data. 

Another possibility is the presence of elevated external 
pressure. In the Galactic center, Oka et al. (1998, 2001) 
find GMCs that have extremely large velocity dispersions 
- approximately a factor of 5 larger at a given radius 
than the clouds in the Galactic disk. Oka et al. (1998) 
argue that the large velocity dispersions from turbulence 
keep the clouds in equilibrium with high external pres- 
sures that arise from hot gas and/or magnetic fields in 
this region. The overlap region of the Antennae may also 
be experiencing higher pressures as the region of con- 
tact between two gas-rich galaxies undergoing a major 
merger. The shocks produced under such violent con- 
ditions, especially in the overlap region, could generate 
high interstellar pressures and high turbulence. Higher 
J transitions of CO are more sensitive to shocked gas, 
and hence the higher velocity dispersions of the Ueda 



et al. clouds may be due to the increased sensitivity of 
CO (3-2) to collisions and turbulence. 

Bolatto et al. (2008) attribute the GMCs that have 
higher velocity dispersions in their sample to blending 
of multiple clouds. This, too, is a distinct possibility in 
our data, given the complex geometry and kinematics of 
the overlap region in the Antennae. Recent work has 
shown that clumpf ind has difficulty accurately measur- 
ing cloud parameters in crowded fields (e.g., Sheth et al. 
2008; Pineda et al. 2009). A few of our spectra in Figure 
5 show hints of multiple peaks that could be interpreted 
as incomplete deblending by clumpf ind. Nevertheless, 
as argued in § 3.5, we do not believe that blending seri- 
ously compromises our results, although it is difficult to 
rule out this effect entirely. 

At this point, it is unclear what is the primary cause 
for the larger velocity dispersions — higher gas densities, 
elevated external pressure, collisions, turbulence, cloud 
blending, or a combination of these factors. Observations 




Fig. 6. — : Molecular gas mass estimated from the CO(2- 
1) intensity versus the virial mass (M v i r = 5ra^/a g G) esti- 
mated from the radius r and velocity dispersion a v found by 
clumpf ind. We assume a geometrical factor a g of 5/3, which 
corresponds to a density profile scaling as r~ 2 . Arrows indi- 
cate upper limits in radius and/or velocity. We show in the 
upper left corner the distance each point would move to the 
left if Xco was 2 and 5 times smaller, and to the right if 
Aco was 2 times larger than our assumed value of Xco = 
3xl0 20 (Kkms -1 ) -1 . Solid line demarcates the one-to-one 
correspondence between the molecular gas and virial mass. 
Most of our points fall on or below the line, suggesting that 
the clouds are gravitationally bound and some may be col- 
lapsing. 

at much higher angular resolution, especially of optically 
thin CO species (C 17 0, C 18 0, 13 CO) and other high- 
density tracers (e.g., HCN, HCO + ) will help distinguish 
between these different possibilities. 

3.7. The Mass Function 

We plot the differential mass function (dN / dM mo \) of 
the clouds found in the SMA+PdBI data in Figure 9, 
with bins of log M mo i/M Q =0.2. The slope of the 
ordinary least squares fit (a) corresponds to a power 
law of the mass function in linear space (dA^/dM mo i oc 
M^ ol ). The improved sensitivity of our observations al- 
lows us to extend the mass function down to lower masses 
than previous work, down to a 5 a completeness limit 
of 3.8x 10 5 M Q 5 . By comparison, Wilson et al. (2003) 
reached a completeness limit more than 10 times higher, 
5x 10 6 M Q (6.7x 10 6 M Q if they adopted a distance of 
22Mpc). 

One of the most remarkable features of the mass 
function is that there is a distinct break around 
log M mo \/MQ w 6.3 — 6.7, which coincides with the 
most massive Galactic clouds (Solomon et al. 1997). This 
break appears to be robust with respect to the detection 
threshold, showing up in both the 2 and 3 <7 c h a n levels. 
Detailed examination of the (u, v) coverage and ampli- 
tude versus (u, v) distance of our SMA and PdBI obscr- 

5 We estimate the 5c completeness limit following the method 
described in Wilson et al. (2003). 
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Fig. 7. — : Relationship between cloud radius (r) and veloc- 
ity dispersion (a v ). Clouds that were unresolved spatially or 
in velocity are assigned upper limits of 80 pc (left-pointing 
arrow) and 4.9 km s - (downward arrow), respectively. Filled 
circles represent our data. Crosses represent GMCs from 
nearby galaxies (NGC 205, NGC 1569, NGC 3077, NGC 4214, 
NGC 4449, NGC 4605, IC 10, SMC, LMC, M33, M31) from 
Bolatto et al. (2008). Open circles represent clouds from 
the Antennae found by Ueda et al. (2012) in CO(3-2). The 
solid line demarcates the typical Milky Way CMC relation- 
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Fig. 8. — : Average H2 number density, roughly esti- 
mated as MH 2 /(§7rr ), as a function of cloud diameter L. 
Crosses represent extragalactic GMCs from Bolatto et al. 
(2008). The solid line denotes the relation (n(H 2 ))(cm~ 3 ) = 
3400 L(pc)" 11 found by Larson (1981). 
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vations does not reveal any gaps that could introduce a 
break in the mass function. The break corresponds to 
the "gap" in the distribution of masses shown in Figure 
3 and is the basis for our motivation to divide the cloud 
population into two subsamplcs. 

For the sake of concreteness, we chose the threshold to 
be log M mo i/M = 6.5. We perform two fits: one in the 
range 4.7 <log M mo i/M < 6.5, and a second in the range 
6.5 <log M mo i/M Q < 8.5. We find a = -1.39 ± 0.10 for 
the first fit and a = —1.44 ± 0.14 for the second fit. 
Figures 3c and 9 indicate that the mass distributions of 
the two populations may be log-normal, dropping off at 
the lower-mass end of each range. If we exclude the two 
end points (log M mo i/M Q = 4.8 and 6.5) where the two 
distributions appears to turn over, the slopes for both 
the low- and high-mass ranges are then steeper, with a = 
-1.52 ± 0.08 and a = -1.52 ± 0.14, respectively. These 
slopes are all statistically consistent with a = — 1.4 ± 0.1 
found by Wilson et al. (2003). 

Note that if we use a distance of 19 Mpc to be com- 
pletely consistent with Wilson et al. (2003), the break re- 
mains apparent at log M mo \/MQ ^6.3 and the respective 
fits have slopes a = — 1.64± 0.16 and a = — 1.40 ± 0.15. 
Additionally, we find that if we run clumpf ind with the 
lower detection thresholds, the slopes of the mass func- 
tion at the low-mass end increases slightly (from —1.39 
to —1.68), while the slope of the high-mass end remains 
relatively constant within the errors. We find that the 
break between the low- and high-mass differential mass 
function also remains constant with different detection 
thresholds; it remains log M mo i/M Q w 6.3 — 6.7 for all 
three thresholds that we tested. 

The slopes of our mass function are well-matched to the 
slope of the mass function for young clusters in the An- 
tennae, a w -2.0 (Zhang & Fall 1999; Whitmore et al. 
2010). This is consistent with the good agreement be- 
tween the GMC mass spectrum and the upper range of 
the SSC mass spectrum in M82 (Keto et al. 2005), sug- 
gesting that individual star clusters may be closely re- 
lated to individual clouds. 

Our fit over the mass range 6.5 <log M mo \/M Q < 8.5 
suggests that some of the largest super-giant molecular 
complexes found by Wilson et al. break up into smaller 
clouds in our high-resolution observations. Importantly, 
many of the large clouds remain large despite the im- 
proved spatial resolution of our observations, which is 
sufficient to resolve the larger structures even at our high- 
est resolution. The largest clouds in our data are much 
larger than GMCs in our Galaxy and typical star-forming 
galaxies (1-100 pc; Solomon et al. 1987; Bolatto et al. 
2008). 

Returning to Figure 3 and Table 3, we see that the 
two populations below and above the break have sta- 
tistically different properties. While the clouds with 
6.5 <log M mo i/M < 8.5 tend to be large (r ;> 100 pc) 

and have high velocity dispersions (a v J> 10 kms -1 ), 
the less massive clouds tend to be largely unresolved in 
both radius (r <J 80 pc) and velocity dispersion (o~ v <J 4.9 

kms -1 ). Also intriguing is the spatial distribution of the 
clouds in Figure 2 — the more massive clouds are asso- 
ciated with peaks of intense star formation (§ 4), while 
the population of smaller, less massive GMCs resides in 
more quiescent areas far from the peaks. 



Slope (4.7-6.5): -1.39+/-0. 10 
Slope (6.5-8.5): -1.44+/-0.14 
Slope: -1.4+/-0.1 (Wilson et al. 03) 




Fig. 9. — : Differential mass function for the clouds extracted 
by clumpf ind from our combined SMA+PdBI data cube. 
Filled black squares show the mass function for clouds found 
above 3<7 c han, and open squares show the mass function for 
clouds found above 2cr c han, showing that the break in the 
mass function is robust at both levels. Unsurprisingly, more 
clouds are found at the lower mass range using the 2 <r c h a n 
lower level, resulting in a steeper mass function for that half 
of the spectrum. Horizontal error bars indicate bin sizes; ver- 
tical error bars indicate \J\N) uncertainty in each bin. The 
vertical dotted line demarcates the 5<j completeness limit of 
our data, 3.8 x 10 s M . We fit two (solid) lines to the 3<r c han 
mass function, one for bins with 4.7<log M mo i/M0 < 6.5 
and a second for bins with 6.5<log M mo i/M0 < 8.5. Both 
fits have slopes that are consistent with that found by Wilson 
et al. (2003) (dashed line). 

In summary, these results suggest a bimodal distribu- 
tion of giant molecular clouds/complexes in the Anten- 
nae: a regular, quiescent population that is similar to 
typical GMCs in nearby galaxies and the Milky Way, 
and a second population that appears much larger, more 
turbulent, and unusually massive. The exceptional prop- 
erties of the largest clouds appear connected to the ex- 
treme conditions associated with the interaction zone of 
the merging galaxies. The dividing line between the two 
populations occurs near M mo \ « 3x 10 6 M Q , which, inter- 
estingly, coincides roughly with the upper end of molec- 
ular cloud masses seen in the Milky Way (Solomon et al. 
1997). 

4. COMPARISON WITH STAR FORMATION TRACERS 

We explore the relationship between dense molecular 
gas and star formation by comparing the CO emission 
with various tracers of recent and on-going star forma- 
tion. 

4.1. Super Star Clusters 

As discussed in the Introduction, SSCs are extremely 
young (< few - 50 Myr) , representing very recent star for- 
mation, and are quite prominent in the Antennae. Figure 
10 shows again the velocity-integrated CO(2-l) emission, 
this time with the locations of various SSCs overplotted. 
These clusters are identified as compact sources in HST 



11 



optical images (Whitmore & Schweizer 1995; Whitmore 
et al. 1999; Zhang et al. 2001). All our CO clouds have 
SSCs nearby. However, the SSCs are identified in the 
optical and are not good tracers of embedded or current 
star formation. 

4.2. Current Star Formation 

We also consider the relationship between current star 
formation and the distribution and intensity of molecu- 
lar gas, velocity gradients, and other signatures of shocks 
and compression in our CO data that may provide in- 
sight into the nature of SSC formation. We focus on 
three star formation tracers with the highest spatial res- 
olution available — Ha emission from the recombination 
of gas ionized by young, massive O and B stars, 8 /im 
emission from polycyclic aromatic hydrocarbons (PAHs) 
heated by star formation, and 4 cm radio continuum from 
either H II regions or young supernova remnants. The 
Ha data were taken with the F658N filter on the HST 
Wide Field Planetary Camera 2 (WFPC2) by Whitmore 
et al. (1999), the PAH maps were acquired using the 
Spitzer Infrared Array Camera (IRAC) by Wang et al. 
(2004), and the 4 cm continuum map was taken by Neff 
& Ulvestad (2000) with the Very Large Array (VLA). 

Wang et al. (2004) estimate that ~5% of the 8 /im 
emission comes from diffuse gas heated by older stars, 
so we subtract this component from the 8 /im image fol- 
lowing Pahre et al. (2004), by scaling the 3.6 and 4.5 /im 
channels to match the colors of M0 III stars. Comparison 
between 4 cm continuum emission and 8 /im PAH emis- 
sion shows very good spatial correlation between the two 
tracers (Figure 11), so we elect to show only the 8 /im 
PAH emission in the next figure. Figure 11 also notes the 
location of peaks in the FIR emission at 70/im (from the 
Herschel Photodetector Array Camera & Spectrometer - 
PACS; Klaas et al. 2010), showing good agreement with 
the 8/im PAH and 4 cm continuum emission within the 
5'.'5 resolution of the FIR image. 

Figure 12 shows the velocity-integrated emission (ze- 
roth moment; top), the intensity- weighted velocity (first 
moment; middle), and the intensity- weighted velocity 
dispersion (second moment; bottom) of our combined 
SMA+PdBI CO(2-l) data. Each panel is overlaid with 
contours of Ha (left) or 8 /im PAH (right) emission to 
show the distribution of star formation. 

As discussed in §3, there appear to be multiple over- 
lapping clouds along the line of sight in the CO-bright 
regions of our data. This means that in a straight first 
moment map, much of the velocity information is lost due 
to averaging. Another way to view the velocity distribu- 
tion of the molecular gas would be to consider position- 
velocity slices of emission at each of the star-forming 
peaks (Figure 13). We choose the position- velocity slices 
to show the largest velocity gradients in the CO emission. 

We discuss the prominent star-forming peaks in the 
8 /im map, labeled A-F in the top right panel of Fig- 
ure 12, and consider the molecular gas distribution and 
kinematics at these peaks. 

Peak A: Star formation in this region is seen both in 
Ha and PAH, both of which appear to be slightly offset 
from the peak CO emission. The higher resolution Ha 
image shows that the star formation breaks up into two 
or three peaks, one of which is almost completely offset 
from the CO. Brandl et al. (2009) estimate a SFR of 



0.74 M Q yr _1 for this whole region based on its infrared 
luminosity. This is also one of the brighter regions in 
terms of integrated CO intensity. The position- velocity 
diagram in Figure 13 shows two or three clouds spatially 
separated by ~3" (~321pc). The spread in gas veloc- 
ity is very large, ranging from 1340 to 1540 kms -1 , but 
there is only the slightest hint of a velocity gradient for 
the upper left CO cloud. Peak A is associated with wa- 
ter maser emission (Brogan et al. 2010), indicating very 
young star formation. 

Peak B: One of the fainter star- forming peaks based on 
the PAH and 4cm emission, this region barely shows up 
in the Ha imaging. The integrated CO emission is also 
one of the weakest, although there may be two clouds 
at around the same velocity, spatially separated by ^3'.'5 
(~375pc). 

Peak C: A region where the Ha emission may be mostly 
obscured, the 8 /im PAH peak is slightly offset from the 
CO peak. This is one of the most highly obscured areas in 
the overlap region, with Ay = 10.3 — 11.8 mag (Snijders 
ct al. 2007; Mengel et al. 2005). The SFR in this region 
is 0.66 Moyr" 1 (Brandl ct al. 2009). The velocity dis- 
persion of the CO is extremely large, >300kms _1 , with 
at least two clouds with large velocity gradients, seen in 
both the middle panel of Figure 12 and the upper right 
panel of Figure 13. Near- infrared spectroscopic imaging 
by Herrera et al. (2011) also find extremely broad line 
widths in the warm H2 emission, which they attribute to 
shocks driven by gas dynamics on large scales. A com- 
pact, warm H2 source observed in this region is thought 
to be associated with a massive cloud on its way to be- 
coming a SSC (Herrera et al. 2011). 

Peak D: This is one of the brightest peaks in 8 /zm 
PAH and in CO(2-l). The large amount of dust and gas 
present explains why this peak is almost completely ob- 
scured in the Ha image. Brandl et al. (2009) estimate the 
SFR for this region to be large, almost 2 M & yr _1 . This 
region also contains a water maser (Brogan et al. 2010). 
Multi-wavelength spectral energy distributions suggest 
that some of the youngest clusters in the Antennae re- 
side here (Zhang et al. 2010). The position- velocity dia- 
gram shows at least two clouds 6 with large velocity gradi- 
ents, possibly in the process of merging. The molecular 
gas spans a large range in velocity, ^200kms _1 , with 
~80kms _1 separation between the two clouds. 

Peak E: This peak has the brightest Ha and second 
brightest 8 /zm PAH emission, but it contains surpris- 
ingly small amounts of CO. Brandl et al. (2009) estimate 
a SFR of 1.8 Mq yr _1 for this region from its infrared lu- 
minosity. Figure 10 shows that the youngest SSCs found 
by Whitmore et al. (1999) are clustered around this re- 
gion. This suggests that peak E is a site of very recent 
and vigorous star formation, but given the lack of molec- 
ular gas present, star formation will not be able to con- 
tinue at the same rate for much longer. The distribution 
of the remaining molecular gas is also interesting; the 
curved shape in the position-velocity diagram is sugges- 
tive of gas infalling toward the center. 

Peak F (NGC 4039): This region contains the nucleus 



6 Note that the contours of CO emission are scaled from the 
faintest cloud, so the lowest contour here is 2x brighter than the 
lowest contour of peak B. 
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Fig. 10. — : Velocity-integrated CO(2-l) emission from the overlap region of the Antennae, with the same contours as Figure 
2. Symbols mark the locations of SSCs of different ages identified by Whitmore & Schweizer (1995) and Whitmore et al. (1999), 
as well as red clusters from Zhang et al. (2001), both of which appear to preferentially surround the CO peaks. Note that the 
red clusters from Zhang et al. (2001) are candidates for the youngest clusters (~5Myr), as the reddening is thought to come 
from dust obscuration. 

Our comparison of the molecular gas distribution and 
kinematics with different star formation tracers (8 /jm 
PAH, 4 cm continuum, and Ha) shows that the molecu- 
lar gas at star-forming peaks have large velocity disper- 
sions, typically spanning velocities in the range of 200- 
300kms _1 . In addition, we find that the molecular gas 
at all but two of the star-forming peaks (B and E) have 
large velocity gradients. 

Keto et al. (2005) find that the clouds with active star 
formation in the starburst galaxy M82 are experiencing 
shock-driven compression, with both lateral and line-of- 
sight velocity gradients (i.e., inverse P Cygni profiles) 
within the individual clouds. These observations suggest 
that individual massive star clusters may have formed 
from individual GMCs crushed by a sudden galactic-scale 
increase in external pressure. While we do not observe 
signatures of inverse P Cygni profiles in our data (the 
presence of multiple clouds along the line of sight may 
blur out such signatures) , the large velocity gradients we 



of NGC 4039, one of the galaxies in the merging An- 
tennae pair. Associated with one of the water masers 
found by Brogan et al. (2010), the age of star clusters 
in this area appears quite young, less than 6.3 Myr (Bas- 
tian et al. 2009). The SFR from the infrared luminosity, 
however, appears much lower than in the other peaks in 
the overlap region (-0.33 M yr _1 ; Brandl et al. 2009). 
Activity from an active galactic nucleus, if any, is ex- 
tremely weak, and the infrared emission appears to be 
dominated by purely young star formation (Brandl et al. 
2009). The distribution of the Ha emission is quite dif- 
ferent and offset from the 8 /im PAH emission. There are 
at least two large clouds in the CO(2-l) emission in this 
region, spanning a velocity range of — 200 km s^ 1 . One of 
the clouds (right) show a fairly obvious velocity gradient, 
while the other (left) does not. The peaks are separated 
spatially by ~4" (-428 pc). 

4.3. Implications 
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Fig. 11. — : Contours of 4cm continuum (black; Neff & Ul- 
vestad 2000) and 8/xm PAH emission (red; Wang et al. 2004) 
in the overlap region of the Antennae, showing very good 
agreement between the two. Contour levels range from 5%- 
95% of peak intensity in steps of 5% for the 4 cm continuum 
emission and 10%— 90% of peak intensity in steps of 10% for 
the 8 fiva. PAH emission. FIR emission peaks from Herschel 
PACS 70/im imaging are marked with crosses, consistent with 
the 8/im PAH and continuum peaks within its 5'.' 5 resolution 
limit. Note that the absence of a FIR counterpart for the 
central 8/im PAH/4 cm continuum peak is due to the coarser 
resolution of the 70/im imaging blending the two peaks to- 
gether. 

observe in some of the star-forming peaks may reflect 
compressive shocks. 

Recent work provides additional observational evi- 
dence supporting the pressure-triggered formation sce- 
nario. As mentioned in the previous section, detailed 
examination of near-infrared emission in Peak C sug- 
gest that the warm H2 emission is excited by large-scale 
shocks (Herrera et al. 2011). Based on the ages of star- 
forming regions, Zhang et al. (2010) propose two paths of 
shock-induced sequential star formation, perhaps propa- 
gated by the collision of GMCs, across the overlap region 
of the Antennae. Both paths start from the center of the 
overlap region (around Peak D in Figure 12) and move 
outwards to the north-west and south-east regions. This 
is well-correlated with the region with the highest star 
formation efficiency (Gao et al. 2001) and FIR emission 
(Klaas et al. 2010). 

Additionally, Whitmore et al. (2010) find evidence for 
shock-triggered sequential cluster formation on a small 
scale, within the individual star-forming regions, based 
on the spatial distribution of clusters of different ages. 
One example of this sequential cluster formation high- 
lighted by Whitmore et al. (2010) involves Peak E in our 
work, where an older cluster (10-50 Myr) to the right 
may have triggered the formation of much younger star 
clusters to the left (Figure 10). And finally, simulations 
of the merger between NGC 4038 and NGC 4039 suggest 
that compressive tides occur frequently and play an im- 
portant role in the evolution of the Antennae (Renaud 
et al. 2009). 



All this clearly ties in with our results from § 3. The 
largest clouds are associated with regions of intense in- 
frared emission, active star formation, and high star for- 
mation efficiencies (Evans et al. 1997; Mirabel et al. 1998; 
Gao et al. 2001; Brandl et al. 2009; Klaas et al. 2010), 
where the molecular gas appear to have large velocity 
dispersions and gradients (§ 4.2). These clouds arc differ- 
ent than typical GMCs in nearby star-forming galaxies. 
While they appear to be gravitationally bound, many of 
these clouds are larger, more massive, and have higher ve- 
locity dispersions than expected from typical GMC scal- 
ing relations. Given that they reside in such an extreme 
environment, some combination of higher gas densities, 
elevated external pressure, turbulence, and cloud colli- 
sions is most likely responsible for the formation of these 
massive clouds, resulting in the bimodal distribution we 
observe in the cloud mass spectrum. 

5. SUMMARY 

In this paper, we present high-resolution SMA and 
PdBI CO(2-l) imaging of the intensely star-forming 
overlap region in the Antennae galaxies in order to study 
the properties of molecular gas that appears to be form- 
ing SSCs. 

Using clumpfind, we resolve the CO(2-l) emission 
into a large number of clouds. This extends the cloud 
mass function down to lower masses than previous work 
by Wilson et al. (2003), down to a ha completeness limit 
of 3.8xl0 5 M Q . A break in the cloud mass function 
around log M mo i/M Q rj 6.5 suggests two distinct popu- 
lation of clouds. Fitting a broken power law to the mass 
function, we find that dN/dM mo i ex M ~1{ for the low- 
mass cloud population and dN/dAf, 



Dl oc 



M~T for 



the high-mass cloud population, close to the M mass 
function found for SSCs in the Antennae. We find that 
the two populations have statistically different proper- 
ties, and appear to populate different areas in the overlap 
region. 

The low-mass population of clouds (log M mo i/M Q < 
6.5), often below our spatial and velocity resolution limit, 
appears to have properties similar to GMCs in typi- 
cal star-forming galaxies. These clouds tend to be dis- 
tributed along the outskirts of the overlap region, away 
from sites of intense star formation. 

The high-mass cloud population (log M mo i/M Q > 
6.5), is associated with sites of peak star formation. 
These clouds appear to be much larger in size (median 
r = 170 pc) and mass (median log M mo i/M Q = 7.2) 
compared to typical GMCs. While we find that these 



clouds appear to be in virial equilibrium (M v 



M mo i), 



they have higher velocity dispersions in the size-linewidth 
relation than expected for typical GMCs. 

Comparison between our CO data and star formation 
traced by 8 /im PAH/4 cm continuum/Ha imaging show 
good correlation between the distribution of molecular 
gas and current star formation. Consistent with obser- 
vational evidence of pressure-triggered star formation in 
the region, position-velocity slices of the CO data re- 
veal large velocity gradients and dispersions in the gas at 
the star-forming peaks, indicative of an active, turbulent 
environment where compressive shocks may dominate. 
These results suggest that factors such as elevated exter- 
nal pressure, increased turbulence, and cloud collisions 
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Fig. 12. — : Color-scale images of the velocity-integrated emission (zeroth moment), intensity-weighted velocity (first moment), and 
intensity-weighted velocity dispersion (second moment) from our combined SMA+PdBI CO(2-l) data from top to bottom, respectively. 
The first and second moments are masked below 10% of peak PAH intensity to suppress the noise in regions of weak CO emission. Color 
bars on the right indicate the range displayed in each panel; left and right panels are identical. The left panels are overlaid with contours of 
Ha emission (Whitmore et al. 1999), and the right panels are overlaid with contours of 8/xm PAH emission (Wang et al. 2004), showing the 
distribution of star formation. Contour levels range from 0.02-1.57 counts s _1 in equal steps of 0.24 dex for the Ha emission and 10%-90% 
of peak intensity in steps of 10% for the 8 /im PAH emission. The star-forming peaks discussed in the text and Figure 13 are labeled A-F 
in the top right panel. 
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Fig. 13. — : Position-velocity diagrams of CO(2-l) emission (upper and lower panels) from our combined SMA+PdBI data, with cuts 
made at the star-forming peaks as traced by 8 fira PAH emission (center panel in heat scale; the cuts are demarcated by white lines and 
labeled A— F to match Figure 12). Cuts are chosen to show the largest velocity gradient in the CO emission. The top and bottom of each 
slice corresponds to the left and right side of each panel, respectively. Contour levels are 20, 35, 50, 65, 80, and 95% of the peak intensity 
of the cut with the weakest emission (Panel C, 0.47 Jy beam -1 ), multiplied by a factor noted in each panel to better display the cuts with 
stronger emission. Color scale of each image is also adjusted accordingly. The position- velocity diagrams show that the structure of the 
molecular gas is often extremely complex at the star-forming peaks, showing multiple peaks and velocities spanning up to 300 kms -1 , 
suggesting that this is a very dynamically active region. 
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are most likely responsible for the apparent high velocity 
dispersions and gas densities, resulting in the formation 
of the high-mass clouds and complexes we observe in the 
overlap region. These high-mass clouds and complexes 
may, in turn, form the SSCs that are so prevalent in the 
Antennae. 

While our observations have allowed us to probe the 
properties of the molecular gas at much higher angu- 
lar resolution and sensitivity than previous work, we are 
still limited to size scales much larger than that of typ- 
ical star clusters (< 10 pc). With the advent of the full 
extended array with ALMA, sensitive and very high- 
resolution (<0'.'l) imaging of the ISM in the Antennae 
will show whether the large, massive clouds we observe 
will begin to break down into smaller objects and steepen 



the cloud mass function. 
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TABLE 2 
Cloud Catalog 



O2000,<52000 ^rad V d v log M mol log M vir 

(kms -1 ) (pc) (kms -1 ) (logM ) (logM©) 
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12 h 01 m 54 s .60 


-18°52'45'.'5 


1601.1 


225 


9.9 


7,4 


7.2 


28 


12 h 01 m 53 s .20 


-18°53'16'.'4 


1679.1 


266 


16.4 


7.5 


7.7 


29 


12 h 01 m 54 s .45 


-18°52'44'.'9 


1620.6 


161 


12.4 


7.2 


7.2 


30 


12 h 01 m 55 s .00 


-18°52'53'.'0 


1562.1 


89 


10.1 


7.0 


6.8 


31 


12 h 01 m 53 s .66 


-18°53'12'/2 


1571.8 


141 


7.7 


7.0 


6.8 


32 


12 h 01 m 53 s .64 


-18°53'll'/9 


1581.6 


119 


11.9 


7.0 


7.1 


33 


12 h 01 m 54 s .89 


-18°52'45'.'5 


1381.6 


118 


9.1 


7.0 


6.8 


34 


12 h 01 m 55 s .65 


-18°52'51'/5 


1537.7 


79 


5.6 


6.7 


6.2 


35 


12 h 01 m 53M6 


-18°53'20'.'9 


1762.0 


113 


9.5 


7.0 


6.9 


36 


12 h 01 m 54 s .55 


-18°52'45'.'2 


1615.7 


94 


12.9 


7.0 


7.0 


37 


12 h 01 m 54 s .91 


-18°53'07'.'l 


1645.0 


216 


10.4 


7.3 


7.2 


38 


12 h 01 m 53 s .92 


-18°53'18'.'2 


1703.5 


211 


7.7 


7.0 


6.9 


39 


12 h 01 m 54 s .40 


-18°52'49'.'7 


1488.9 


55 


10.7 


6.9 


6.6 


40 


12 h 01 m 54 s .87 


-18°52'59'.'6 


1610.9 


204 


14.4 


7.2 


7.5 


41 


12 h 01 m 54 s .81 


-18°52'47'.'6 


1386.5 


261 


18.9 


7,4 


7.8 


42 


12 h 01 m 54Ml 


-18°53'26'.'3 


1606.0 


53 


4.2 


6.4 


5.8 


43 


12 h 01 m 55 s .23 


-18°53'09'.'5 


1630.4 


168 


8.7 


6.8 


6.9 


44 


12 h 01 m 54 s .00 


-18°53'09'.'5 


1498.7 


70 


7.3 


6.7 


6.4 


45 


12 h 01 m 54 s .57 


-18°52'46'.'4 


1430.4 


229 


11.7 


7.2 


7.3 


46 


12 h 01 m 55 s .25 


-18°53'09'.'2 


1640.1 


75 


7.2 


6.6 


6,4 


47 


12 h 01 m 54 s .74 


-18°52'40'.'4 


1513.3 


<80 


12.1 


6.6 


<6.9 


48 


12 h 01 m 54 s .28 


-18°52'50'.'3 


1532.8 


134 


11.9 


6.9 


7.1 


49 


12 h 01 m 54 s .55 


-18°52'45'.'8 


1396.3 


71 


5.9 


6.5 


6.2 


50 


12 h 01 m 55 s .33 


-18°52'43'.'7 


1645.0 


176 


12.0 


6.8 


7.2 


51 


12 h 01 m 54 s .76 


-18°52'53'.'3 


1527.9 


102 


16.0 


6.8 


7.3 


52 


12 h 01 m 55 s .48 


-18°53'07'.'4 


1601.1 


149 


10.0 


6.8 


7.0 


53 


12 h 01 m 53 s .66 


-18°53'12'/8 


1552.3 


97 


7.0 


6.5 


6.5 


54 


12 h 01 m 55 s .02 


-18°52'52'.'7 


1552.3 


227 


18.1 


7.1 


7.7 


55 


12 h 01 m 54M7 


-18°52'47'.'3 


1513.3 


109 


7.1 


6.6 


6.6 


56 


12 h 01 m 55 s .25 


-18°52'43'.'l 


1537.7 


114 


10.0 


6.6 


6.9 


57 


12 h 01 m 54 s .26 


-18°52'49'.'4 


1523.1 


127 


12.8 


6.8 


7.2 


58 


12 h 01 m 54M7 


-18°52'45'.'8 


1615.7 


<80 


2.1 


6.3 


<5.4 


59 


12 h 01 m 54 s .49 


-18°52'45'.'2 


1406.0 


128 


10.6 


6.6 


7.0 


60 


12 h 01 m 53 s .64 


-18°53'08'.'9 


1527.9 


238 


11.8 


7.1 


7,4 


61 


12 h 01 m 54M9 


-18°52'44'.'6 


1610.9 


<80 


9.4 


6.3 


<6.7 


62 


12 h 01 m 53 s .94 


-18°53'17'/9 


1684.0 


220 


6.2 


6.8 


6.8 


63 


12 h 01 m 54 s .66 


-18°52'4T/3 


1518.2 


35 


7.4 


6.6 


6.1 


64 


12 h 01 m 53 s .28 


-18°53'2T/8 


1737.7 


<80 


11.2 


6.3 


<6.8 


65 


12 h 01 m 55 s .65 


-18°52'41'/9 


1606.0 


144 


11.0 


6.5 


7.1 


66 


12 h 01 m 54 s .38 


-18°52'47'/3 


1527.9 


163 


13.9 


6.9 


7.3 
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TAI 


3LE 2 - 


Continued 










O2000 


feooo 


*>rad 


r 


<T V 


log M mol 


log M vir 








(kms -1 ) 


(pc) (kms -1 ) 


(logM ) 


(logM ) 


67 


12 h 01 m 53 s .24 


-18°53'17'/0 


1762.0 


189 


10.4 


6.8 


7.2 


68 


12 h 01 m 53 s .69 


-18°53'13'/4 


1645.0 


166 


24.0 


6.8 


7.8 


69 


12 h 01 m 54M7 


-18°52'45'.'8 


1601.1 


<80 


6.3 


6.1 


<6.3 


70 


12 h 01 m 54 s .72 


-18°52'53'/3 


1547.5 


99 


10.7 


6.7 


6.9 


71 


12 h 01 m 54M3 


-18°52'47'/3 


1459.7 


<80 


<4.9 


5.1 


<6.1 


72 


12 h 01 m 54M7 


-18°53'17'/0 


1532.8 


1.30 


9.0 


6.5 


6.9 


73 


12 h 01 m 55 s .50 


-18°52'52'.'4 


1474.3 


<80 


1.6 


5.9 


<5.1 


71 


12 h 01 m 53 s .07 


-18°53'18'/2 


1445.0 


<80 


<4.9 


5.6 


<6.1 


75 


12 h 01 m 55 s .59 


-18°52'53'/3 


1493.8 


<80 


2.0 


5.9 


<5.3 


76 


12 h 01 m 55 s .50 


-18°52'52'/7 


1440.2 


<80 


3.5 


5.9 


<5.8 


77 


12 h 01 m 53 s .73 


- 18° 53' 13'.' 1 


1664.5 


<80 


2.0 


5.8 


<5.4 


78 


12 h 01 m 54 s .21 


-18°52'43'.'7 


1415.8 


<80 


<4.9 


5.2 


<6.1 


79 


12 h 01 m 53 s .60 


-18°53'17'/6 


1430.4 


<80 


<4.9 


5.5 


<6.1 


80 


12 h 01 m 54 s .05 


-18°53'18'.'2 


1498.7 


<80 


4.2 


5.9 


<6.0 


81 


12 h 01 m 55 s .38 


-18°53'05'/6 


1567.0 


<80 


2.5 


5.9 


<5.5 


82 


12 h 01 m 54 s .30 


-18°52'41'.'9 


1518.2 


<80 


1.2 


5.7 


<4.9 


83 


12 h 01 m 53M6 


-18°53'09'/5 


1479.2 


105 


<4.9 


6.1 


<6.2 


81 


12 h 01 m 54 s .26 


-18°52'54'/2 


1659.6 


<80 


<4.9 


5.1 


<6.1 


85 


12 h 01 m 54 s .40 


-18°52'44'.'3 


1376.8 


<80 


<4.9 


5.8 


<6.1 


86 


12 h 01 m 55 s .61 


-18°52'51'.'5 


1508.4 


<80 


<4.9 


5.4 


<6.1 


87 


12 h 01 m 53 s .18 


-18°53'15'.'5 


1396.3 


<80 


<4.9 


5.3 


<6.1 


88 


12 h 01 m 54 s .53 


-18°53'24'.'8 


1513.3 


11 


2.1 


6.0 


5.2 


89 


12 h 01 m 55 s .23 


-18°52'41'.'0 


1615.7 


<80 


<4.9 


5.2 


<6.1 


90 


12 h 01 m 53 s .69 


-18°53'10'/7 


1488.9 


<80 


<4.9 


5.5 


<6.1 


91 


12 h 01 m 54 s .24 


-18°52'45'/2 


1367.0 


<80 


<4.9 


5.3 


<6.1 


92 


12 h 01 m 54 s .53 


-18°53'17'/3 


1674.3 


<80 


1.8 


5.7 


<5.2 


93 


12 h 01 m 54Ml 


-18°53'26'/6 


1615.7 


<80 


<4.9 


5.4 


<6.1 


91 


12 h 01 m 53 s .88 


-18°53'08'.'3 


1479.2 


<80 


0.8 


5.6 


<4.5 


95 


12 h 01 m 54 s .85 


- 18° 53' 16'.' 1 


1664.5 


<80 


<4.9 


5.5 


<6.1 


96 


12 h 01 m 53 s .71 


-18°53'13'.'7 


1527.9 


<80 


<4.9 


5.2 


<6.1 


97 


12 h 01 m 53 s .20 


-18°53'08'.'3 


1464.5 


<80 


<4.9 


5.6 


<6.1 


98 


12 h 01 m 53 s .58 


- 18° 53' 10'.' 1 


1732.8 


<80 


<4.9 


5.1 


<6.1 


99 


12 h 01 m 54 s .ll 


-18°52'49'.'l 


1606.0 


<80 


<4.9 


5.3 


<6.1 


100 


12 h 01 m 54 s .ll 


-18°53'26'.'3 


1649.9 


<80 


1.0 


5.4 


<4.7 


101 


12 h 01 m 54 s .34 


-18°52'47'.'9 


1567.0 


<80 


<4.9 


5.6 


<6.1 


102 


12 h 01 m 54M9 


-18°53'18'.'5 


1503.6 


<80 


<4.9 


5.3 


<6.1 


103 


12 h 01 m 54 s .91 


-18°52'38'.'9 


1328.0 


<80 


<4.9 


5.1 


<6.1 


101 


12 h 01 m 53 s .79 


-18°53'20'.'6 


1484.1 


<80 


4.5 


5.7 


<6.0 


105 


12 h 01 m 53 s .73 


-18°53'19'.'7 


1727.9 


<80 


<4.9 


5.1 


<6.1 


106 


12 h 01 m 54M3 


-18°53'26'.'9 


1635.2 


<80 


<4.9 


4.9 


<6.1 


107 


12 h 01 m 54 s .28 


-18°52'51'.'5 


1669.4 


<80 


<4.9 


5.3 


<6.1 


108 


12 h 01 m 54 s .21 


-18°53'll'/6 


1537.7 


<80 


<4.9 


5.2 


<6.1 


109 


12 h 01 m 54M7 


-18°52'49'.'7 


1479.2 


<80 


<4.9 


5.3 


<6.1 


110 


12 h 01 m 55 s .44 


-18°52'45'.'8 


1645.0 


<80 


<4.9 


5.1 


<6.1 


111 


12 h 01 m 53 s .22 


-18°53'12'.'8 


1488.9 


<80 


<4.9 


5.0 


<6.1 


112 


12 h 01 m 54 s .02 


-18°53'16'.'4 


1742.5 


<80 


<4.9 


5.5 


<6.1 


113 


12 h 01 m 54 s .21 


-18°53'21'/8 


1732.8 


<80 


<4.9 


5.2 


<6.1 


111 


12 h 01 m 54 s .28 


-18°52'43'.'4 


1537.7 


<80 


<4.9 


5.1 


<6.1 


115 


12 h 01 m 54 s .72 


-18°53'05'.'9 


1654.8 


<80 


<4.9 


5.5 


<6.1 


116 


12 h 01 m 53 s .18 


-18°53'09'.'2 


1493.8 


<80 


1.3 


5.3 


<4.9 


117 


12 h 01 m 54 s .62 


-18°53'08'.'0 


1381.6 


<80 


<4.9 


5.0 


<6.1 


118 


12 h 01 m 53 s .88 


-18°53'll'/6 


1527.9 


<80 


<4.9 


5.0 


<6.1 


119 


12 h 01 m 54 s .26 


-18°52'44'.'0 


1547.5 


<80 


1.0 


4.9 


<4.8 


120 


12 h 01 m 54 s .98 


-18°52'44'.'3 


1645.0 


<80 


<4.9 


5.3 


<6.1 


121 


12 h 01 m 53 s .43 


-18°53'18'.'8 


1391.4 


<80 


<4.9 


1.9 


<6.1 


122 


12 h 01 m 53 s .69 


-18°53'17'/6 


1718.2 


<80 


<4.9 


4.9 


<6.1 


123 


12 h 01 m 55 s .48 


-18°52'45'.'5 


1620.6 


<80 


<4.9 


5.1 


<6.1 


121 


12 h 01 m 54M9 


-18°52'47'.'6 


1674.3 


<80 


<4.9 


4.9 


<6.1 


125 


12 h 01 m 53 s .54 


-18°53'18'.'2 


1786.4 


<80 


<4.9 


4.9 


<6.1 


126 


12 h 01 m 55 s .61 


-18°52'43'.'7 


1586.5 


<80 


<4.9 


1.9 


<6.1 


127 


12 h 01 m 54Ml 


-18°53'26'.'9 


1537.7 


<80 


<4.9 


1.9 


<6.1 


128 


12 h 01 m 54 s .28 


- 18° 53' 16'.' 1 


1562.1 


<80 


<4.9 


4.9 


<6.1 


129 


12 h 01 m 53 s .92 


-18°53'06'.'5 


1513.3 


<80 


<4.9 


4.9 


<6.1 


130 


12 h 01 m 54 s .09 


-18°53'25'.'l 


1635.2 


<80 


<4.9 


5.1 


<6.1 


131 


12 h 01 m 54M5 


-18°53'10'.'7 


1640.1 


<80 


<4.9 


5.0 


<6.1 


132 


12 h 01 m 54 s .74 


-18°53'05'.'9 


1674.3 


<80 


1.3 


4.9 


<4.9 
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TABLE 2 — Continued 



"2000,<52000 *>rad 1~ °~v log M mol log M vir 

(kms" 1 ) (pc) (fans" 1 ) (logM©) (logM ) 

Note. — Basic cloud properties from clumpfind. Clouds that were unresolved spatially or in velocity arc assigned upper limits of SO pc and 
4.9kms~ , respectively. Note that the velocity is defined using the radio definition (v — cz/(l + z)). 



